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We discuss the effects of cosmic phase transition on the spectrum of primordial gravitational 
waves generated during inflation. The energy density of the scalar condensation responsible for the 
phase transition may become sizable at the epoch of phase transition, which significantly affects the 
evolution of the universe. As a result, the amplitudes of the gravitational waves at high frequency 
modes are suppressed. Thus the gravitational wave spectrum can be a probe of phase transition in 
the early universe. 



Spontaneous symmetry breaking (SSB) often plays 
very important role in high energy physics. For the 
construction of the standard model of particle physics, 
which is currently the most successful model of high en- 
ergy phenomena, the SSB of SU{2)l x U{1)y -> ?7(l)om 
due to the Higgs mechanism (i.e., the electroweak sym- 
metry breaking) is crucial. In addition, at the scale of 
the quantum chromodynamics (QCD), chiral symmetry 
breaking also occurs. If we consider various models of 
physics beyond the standard model, SSBs may occur at 
higher energy scales. For example, if the strong CP prob- 
lem is solved by the Peccei-Quinn (PQ) mechanism 
the PQ symmetry breaking should happen at the PQ 
scale. In grand unified theories (GUTs) [2], the sym- 
metry breaking of Gqut ^ SU{i)c x SU{2)l x U{1)y 
occurs at the GUT scale (where Ggut is the GUT gauge 
group). In supersymmetric models, the supersymmetry 
breaking terms are expected to arise due to the sponta- 
neous breaking of supersymmetry. 

The SSBs in the framework of the standard model, 
i.e., the electroweak symmetry breaking and the chiral 
symmetry breaking in QCD, may be well understood in 
the future by experimental data (in particular, by the 
LHC result), lattice simulation, and so on. However, it 
is difficult to study the SSBs in models beyond the stan- 
dard model because their energy scales are too high to 
be reached by collider experiments. Thus, the physics 
related to those SSBs at high energy scales are hardly 
probed by the existing methods. 

If we consider cosmology, there should exist in the past 
a period of cosmic phase transition related to the SSB. 
In particular, in a period around the phase transition, 
the expansion of the universe may be significantly af- 
fected by the energy density of the fields which cause 
the SSB, which results in a significant deviation from 
the radiation-dominated universe. In the following, we 
will show that information of such an early universe may 
be imprinted in the spectrum of primordial gravitational 
waves (GWs) generated during inflation. This is because 
the GW spectrum is sensitive to the expansion history of 
the universe Importantly, the GW spectrum r riay 

be studied at future experiments such as DECIGO |ig | 
and/or BBO 

In this letter, we discuss the possibility of studying 
cosmic phase transition in the early universe by using 
GW. We will show that the evolution of the amplitudes 



of GW may be significantly affected if phase transition 
happened in the early epoch and hence the information 
on the SSB may be extracted from the spectrum of GW. 

Let us first discuss how the amplitudes of GWs evolve 
in the expanding universe. The tensor perturbation of 
the metric, which corresponds to the degrees of freedom 
of the GW, is defined as 

ds^ = -dt^ + a^{t) {S.,j + 2h,j) dx^dx\ (1) 

where ds is the line element, and the indices i and j run 
1 — 3. The tensor perturbation hij satisfies transverse 
and traceless conditions, h\ = hif-' — 0. Thus, there 
are two physical degrees of freedom (for a fixed value 
of the momentum) , which we denote by -I- and x . In 
the following discussion, it is convenient to work in the 
momentum space, so we define the Fourier amplitude of 
h'ij as 

where Mpi ~ 2.4 x 10^^ GeV is the reduced Planck 
scale, and e"^^ is the polarization tensor which satisfies 

e-^^e|^ = 5\\'. Then, /li^^ satisfies 

l^^^Hi^ + ^-O, (3) 

where the "dot" denotes the derivative with respect to 
time, and k = \k\. Here, we neglect the anisotropic stress 
which is irrelevant for the present study. Although we 
will numerically follow the evolution of the GW ampli- 
tudes, it is instructive to shortly discuss the qualitative 
behavior of the solution of Eq. (jS]). When k <^ aH (out 
of horizon), the last term of the left-hand side is irrele- 
vant and hence h-^ stays (almost) constant. On the con- 
trary, once the mode enters the horizon (i.e., k ^ aH), 

is under oscillation; in this case, (/i|)osc and (/i|)osc 
are approximately proportional to and a"**, respec- 
tively, where (• • • )osc denotes the average for the time 
scale much longer than the oscillation time (but shorter 
than cosmic time). 

During infiation, the quantum fluctuation of the GW 
determines the initial value of h'^-^K The present GW 
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spectrum per log wavenumber interval for fcEw ^ <C 
^RH (where /cew and /crh are comoving wavenumbers 
of modes which enter the horizon at the time of elec- 
troweak phase transition and reheating after inflation, 
respectively) is given by 



n'§^\k)^1.7xl0-'^ro.il ■■ fcEW<fc«fcRH. (4) 
Here, ro.i is the tensor-to-scalar ratio in units of 0.1 and 
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where Tin(fc) denotes the temperature at which the mode 
k enters the horizon, and g^s denote the effective 
number of relativistic degrees of freedom for the en- 
ergy density and the entropy density, respectively, with 
subscript being for the present value. In addition, 
ko = 0.002Mpc~"'^ is the pivot scale and rif is the tensor 
spectral index, which is given by rit = — r/8 in standard 
inflation models. Thus the scale dependence of the pri- 
mordial GW spectrum is very weak as long as r is small 
enough. Hereafter, we neglect the scale dependence for 
simplicity. 

Now let us discuss how the universe expands in the pe- 
riod of cosmic phase transition. During the cosmic phase 
transition, the expectation value of the order parame- 
ter changes from zero to a finite value due to thermal 
effects. The detail of the phase transition depends on 
the physics in the SSB sector. In the present study, we 
model the SSB sector simply by introducing a scalar field 
(j) which plays the role of the order parameter. With 
and X being (real) scalar fields, we consider the following 
scalar potential. 



F(0) 
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where g and h are coupling constants while is the vac- 
uum expectation value of cj)."^^ Here, represents the 
scalar field responsible for the cosmic phase transition; 
(/) = (0 = v^) corresponds to the symmetric (broken) 
phase, while x represents the degrees of freedom in ther- 
mal bath (with temperature T). 

With X being in thermal bath, the free energy of ac- 
quires a term proportional to T^(j)^ around cj) — 0. Thus, 
in the early universe with high enough temperature, the 
symmetry is expected to be restored. As the temperature 
decreases, the negative mass-squared term in the (zero- 
temperature) scalar potential wins over the thermal mass 
term, and the SSB takes place. Here, we expect that the 
evolution of (f> is well governed by the following equation. 



(7) 



where H is the expansion rate of the universe, is 
the decay rate of (j>, and the "dash" denotes deriva- 
tive with respect to (f). In addition, Vt is the poten- 
tial of (f> in thermal bath, which we evaluate as Vr = 
- vlf + ^h{x^)T(b^ where {0}t = tr[Op] (with 
p being density matrix) is thermal average of the oper- 
ator O. We consider the case where the time scale of 
the change of </> is much smaller than that of cosmic ex- 
pansion. Thus, we take p — e^-^^^"^ , where is the 
Hamiltonian in x-sector to calculate (x^)t- 

We evaluate the thermal average by approximating x 
as a free scalar field with the mass squared of hcj)'^ and 
obtain 



^t(0) 



15( 
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(8) 



where ^ = y/p'^ + hcj)^.'^'^ For the following discussion, 
it is instructive to expand the potential around = to 
find 



^T(0) = ^Mr'-TcV 



where 




(9) 



(10) 



Thus, the curvature at = changes its sign at T = Tc. 

Because we are interested in the case where F^ ^ i?, 
the oscillation of (j) decays away with the time scale much 
faster than the cosmic expansion. Thus, we can approxi- 
mate that (j) follows the temporal minimum of Vt- Even 
so, an accurate understanding of the evolution of is 
not straightforward because, in some period, Vt has two 
minima. In particular, at the temperature just above Tc, 
= is not the absolute minimum of the potential, and 
hence the phase transition may proceed with a tunnel- 
ing from (/) = to the absolute minimum [T^j. Whether 
the phase transition completes via the first or second or- 
der phase transition depends on the parameters in the 
model. For ft, ^ 1, the second order phase transition 
precedes the first order one. For h ^ 0(1), the first or- 
der phase transition may take place at the temperature 
T (a few) X Tc. In either case, we simply approximate 
that the phase transition occurs at the time when the 
cosmic temperature becomes Tc for the first time. In ad- 
dition, because the position of the true minimum of the 



Problematic domain wall formations are avoided if one regards 
as a radial part (absolute value) of the complex scalar field, which 
triggers the SSB of a continuous symmetry, such as U(l). Cosmic 
strings in association with the SSB of U(l) are not harmful for 
< 10^5 GeV. 



One may estimate the expectation value of (j> by minimizing the 
free energy. In the present analysis, we approximate that the 
phase transition occurs when the curvature of the potential at 
if> becomes zero. With such an approximation, results of free- 
energy and our procedures agree because the thermal mass of 
the scalar field 4> obtained in two procedures are the same. 
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FIG. 1: Evolution of the product tH as a function of time 
(normalized by t^^^^^) for g = 10~^ (red-solid), 10^'' (green- 
dashed), and 10~* (blue-dotted). Here we take h — \. Notice 
that the figure is independent of v^. 
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FIG. 2: nGw(fc)/n|'^(fc) as a function of k (normalized by 
fcpT) for g = 10"^ (red-solid), 10"^ ( green-dashed), and 10 
(blue-dotted). Here we take h — 1. Notice that the figure is 
independent of 



potential at T = well agrees with that at T = 0, we 
approximate the energy density of the (j) sector as 



P<f> 



V{0) 




t < tpT 
t > tpT 



(11) 



where tpT is the time of the phase transition. 

With the above approximation, the evolution of the 
scale factor a is governed by 



Piad + P<p 



(12) 



where prad is the energy density of the radiation compo- 
nent; it evolves as 



Prad + 4iJprad = V{0)S{t - tp^) , 



(13) 



and is related to the cosmic temperature as prad = 

2 . 

^g*T , with being the effective number of massless 
degrees of freedom. In our analysis, we take the standard- 
model value of g*, which is 106.75. Although x may also 
contribute to 17*, such a contribution is so small that we 
can safely neglect it. 

In Fig. [1] we plot the product of the time t and the 
expansion rate H as a function of t for several values 



(9=1) 
PT 



which is the 



of g. (In the plot, t is normalized by t 
time of the phase transition for the case of 5 = 1, to make 
the figure independent of v^.) The product tH is equal 
to ^ if the universe is dominated by radiation. We can 
see that the evolution of the universe at i ^ ipT deviates 
from that of radiation-dominated universe as g becomes 
smaller. This behavior can be easily understood from the 
relation [p0/piad]t=tpT ~ 0{h? / g^g). For smaller the 
potential energy of (j) at the origin tends to dominate the 



universe before the phase transition. In the small g limit, 
a brief period of inflation takes place . 

Once the evolution of the scale factor is understood, we 
can easily solve Eq. ^ to obtain the present spectrum of 
the GW. In Fig. [21 we plot the present GW spectrum as 
a function of k. Assuming that Tc is much higher than 
the electroweak scale, we normalize JIqw as 



f^Gw(fc«fcpT) = r![f^^(fc) 



(14) 



where kp^p = apx^^PT is the wavenumber of the mode 
which enters the horizon at the time of the phase transi- 
tion. 

One can see that the GW spectrum with k > fcpx is 
suppressed. With the present approximation, the follow- 
ing relation holds. 
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GW 
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Prad(Tc)+F(0)^ 



(15) 



where R is the reduction rate of the high-frequency GW 
spectrum due to the phase transition. The right-hand 
side of Eq. ([T5|) depends only on the combination of g/fi^, 
and is independent of w^. For g/h"^ = Ix lO^'* (3 x lO^'*, 
1 X 10-3, 3 X 10-3, 1 X 10-2), R is given by 0.25 (0.50, 
0.76, 0.91, 0.97). If a short period of inflation occurs with 
sufficiently small g/h'^, the spectrum of GWs which enter 
the horizon during such a period is proportional to k~*. 

In order to discuss the possibility of studying the cos- 
mic phase transition using GWs, it is necessary to under- 
stand the present frequency of the mode with k ^ kp^. 
(The comoving wavenumber is related to the present fre- 
quency as / = fc/27rao, with oq being the present scale 
factor.) Let us define 



PT 



5g 



1/4 



(16) 
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which is the temperature just after the phase transition. 
Then, the present frequency of the mode with k — fcpT 
is given by 

- ^ ( 10^) ■ ^''^ 

We have seen that g/h^ -C 1 when 1 — R becomes sizable. 
In such a case, we obtain 

[/pt],/..«i - 0.50 Hz X (^^) . (18) 

Finahy we discuss the possibihty for detecting charac- 
teristic features of the phase transition in the GW spec- 
trum. For this purpose, we approximate the sensitivity of 
the future space interferometers such as DECIGO/BBO 
with correlation analysis of 1 year [HI, [3 [3 as 

r nt^^\h/f) : f<h, 

ac„s(/) = <^ n^^:s\f/hf : h<I<h. 

[ ni^s\j2/hnf/f2f : /2</, 

(19) 

with niTns^ = 3 X 10-l^ /i = 0.3 Hz, and /2 = 2 Hz. 
Then, we expect that the modulation in the GW spec- 
trum due to the phase transition is in the range of detec- 
tor sensitivity if (i) (1 - R)^GwifpT) > ^^sons(/pT), and 
(ii) R^lGwifpT) > risons(max{/pT, /i}). Condition (i) 
ensures that the drop-off of ficw is larger than the sen- 
sitivity, while condition (ii) means that the GWs with 



k ^ kpT are observable. For (r, /pt) = (0.1,0.1 Hz), 
(0.1,1 Hz), and (0.01,0.1 Hz), for example, the condi- 
tions (i) and (ii) are satisfied when 0.005 < R < 0.98, 
0.17 < R < 0.83, and 0.05 < R < 0.86, respec- 
tively. Broader regions will be explored in the ultimate- 
DECIGO [13, [3| , where sensitivities will be improved by 
orders of magnitude. Notice that, because of the stochas- 
tic background from white dwarf binaries, it will be dif- 
ficult to extract the signal of cosmic phase transition in 
the GW spectrum if /pt < 0.1 Hz [i3]. 

As a final remark, the scalar field dynamics associated 
with phase transitions produces GWs of flat spectrum 
[17H19l |. This contribution is small enough to be ne- 
glected for the intermediate scale phase transition with 

^ 10® GeV, which we are interested in (see Eq. ([TSl) ). 

In summary, we have argued that the spectrum of GW 
can be a useful probe of the cosmic phase transition in the 
early universe. So, if GWs with sizable amplitudes (i.e., 
the sizable value of the tensor-to-scalar ratio parameter 
r) are observed by the measurement of -B-mode polar- 
ization of the cosmic microwave background, we have a 
good chance of studying the cosmic phase transition in 
the early universe with precise observations of primordial 
GWs. 
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